Minimal seesaw model with S4 flavor symmetry 
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We discuss a neutrino mass model based on the ^4 flavor symmetry within the minimal seesaw 
framework, in which only two right-handed neutrinos are introduced and transform as 2 under 
5*4. Although the model contains less free parameters compared to the typical seesaw models, it 
provides a successful description of the observed neutrino parameters, and in particular, a nearly 
tri-bimaximal mixing pattern can be naturally accommodated. In addition, the heavy right-handed 
neutrino masses are found to be non-degenerate, while only the normal hierarchical mass spectrum 
is compatible with experiments for Ught neutrinos. 



I. INTRODUCTION 



In view of the compelling experimental evidence on 
neutrino oscillations, the origin of neutrino masses and 
Icpton flavor mixing emerges as one of the most funda- 
mental issues in particle physics. Since neutrinos are 
massless particles in the standard model (SM) of par- 
ticle physics, a broad class of models extended the SM 
have been proposed in order to accommodate massive 
neutrinos. The seesaw mechanism [H-Q turns out to 
be among the most attractive extensions of the SM in 
virtue of its natural explanation of tiny neutrino masses. 
In the canonical type-I seesaw model, three heavy right- 
handed neutrinos are introduced besides the SM particle 
contents, while a Majorana mass term Mji is assumed, 
which is not subjected to the scale of electrowcak sym- 
metry breaking scale, i.e., Aew ^ 100 GcV. The light 
neutrino mass scale is then strongly suppressed with re- 
spect to Aew due to the heavy right-handed neutrino 
masses. 

In general, the type-I seesaw model is pestered with too 
many model parameters, and therefore, fails to predict 
the lepton flavor mixing pattern as well as the light neu- 
trino mass spectrum. For example, in case of the simplest 
type-I seesaw model with three right-handed neutrinos, 
there are in total fifteen free parameters in the Dirac mass 
matrix together with three unknown mass eigenvalues of 
heavy Majorana neutrinos, whereas the light neutrino 
mass matrix contains only nine physical parameters, in- 
dicating a lack of valuable predictions. Note that, in the 
most economical type-I seesaw model, i.e., the minimal 
seesaw model (MSM) [ol-fl^. one could introduce only 
two right-handed neutrinos, whereas the observed neu- 
trino mass hierarchy and lepton fiavor mixing can be well 
interpreted. Such a minimal extension of the SM greatly 
reduces the number of free parameters, and hence is very 
predictive. For instance, in the MSM, one of the light 
neutrinos should be massless since is of rank 2, which 
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indicates that m; ~ 0.05 cV in the normal hierarchy 
case while rrii ~ 0.1 eV in the inverted hierarchy case 
with rrii being the light neutrino masses. In case that 
future cosmological observations set more stringent con- 
straints on the summation of the light neutrino masses, 
the MSM would then be the most plausible underlying 
model. 



Recently, plenty of models extended the gauge group 
with flavor symmetries are studied in order to understand 
the lepton flavor mixing. In particular, the experimen- 
tally favored tri-bimaximal mixing pattern jl3l4l5| can 
be naturally realized in many flavor symmetry models. 
It is therefore interesting to investigate if the neutrino 
masses and mixing can be realized in the MSM based on 
certain flavor symmetries. Now that there are only two 
right-handed neutrinos in the MSM, the symmetry group 
G f should contain at least one two-dimensional represen- 
tation, if two right-handed neutrinos are located in the 
same multiplet of G/. In addition, a three-dimensional 
representation should be employed in order to accommo- 
date three generations of charged leptons in a natural 
way. In this sense, the permutation group appears as 
an attractive candidate for the MSM, since it is one of the 
smallest discrete groups containing one-, two- and three- 
dimensional representations. Similar models of the 5*4 
flavor symmetry within the canonical seesaw framework 
have been intensively studied in the literature [l6l - l46j . 

In this work, we consider the MSM based on the 5*4 
flavor symmetry. In particular, we shall show that our 
scheme is rather compact whereas it is compatible with 
the experimental observation, i.e., the tri-bimaximal mix- 
ing pattern could be easily accommodated. The remain- 
ing parts of the work is organized as follows: In Sec. [ill 
we present the main content of our model, and formu- 
late the general expressions of the lepton mass matrices. 
One interesting example is given in order to show how 
the tri-bimaximal mixing is realized. The information 
on the Higgs potential is also briefly discussed. Then, in 
Sec. llin we perform a detailed numerical analysis, and il- 
lustrate the main results obtained in the model. Finally, 
our conclusions are presented in Sec. IIVI 
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II. THE MODEL 

The discrete group S4 is the permutation group of four 
distinct objects, which contains 24 group elements and 5 
irreducible real representations. Among the five repre- 
sentations, two are one-dimensional (li and I2), one is 
two-dimensional (2), and two are three dimensional (3i 
and 32). The group properties, i.e., the Kroneckcr prod- 
ucts and the Clebsch Gordan coefficients, can be found 
in the appendices of Ref. [l^l • 

The total symmetry of our model is then chosen to be 

G = SU{3)c(E>SU{2)(g)U{l)Y(E>S4, (1) 
under which the lepton content in our model is placed as 



(1,2,-1)(32), (2) 
^ (l,l,-2)(32), (3) 
(1,1,0)(2). (4) 

the minimal seesaw model is con- 
sidered whereas the two right-handed neutrinos are as- 
signed to the trivial representation of ^4. Furthermore, 
the Higgs assignments in our model are given by 

(l,2,-l)(li), (5) 



L 

Note that in Ref. Ei 



('/'I, 02) 

(6,6,6) 

(Xl,X2) 



(1,2,-1)(2), 

(l,2,-l)(3i), 

(1,1,0)(2), 



(6) 
(7) 
(8) 



where the SUVi) doublets Higgs fields arc in analogy to 
these in Ref. [201, whereas an additional SU{2) singlet 
Higgs X is introduced. We will show later on that x is 
crucial to ensure the corrected prediction on the neutrino 
mixing angles as well as the light neutrino masses. Note 
that we mainly focus our attention on the lepton flavor 
mixing, and hence do not include the quark sector in our 
discussions. A simple way to contain the quark mixing 
in our model is to make a naive assumption that all the 
quarks belong to the identity representation, i.e., li, and 
then the quark flavor mixing and masses can be obtained 
via the standard Yukawa couplings to 0o- 

By using the group algebra of 5*4, we can write the 
invariant Yukawa couplings for leptons as 

C = ao {Licr + L2fJ-R + Lstr) (j)o 

+ ai VS {L2fJ.R - LsTr) (t>i 

+ {-2LieH + L2fJ-R + L3Tfl)02] 
+ a2 [{L2TR + LjfiR) ^1 + {LiTr 

LiVRid + \^-L2VRi 

iL^VR2 6 



^R2^RV X2J 
(9) 




Lscr) 6 



3L2l^R2 I 6 



J^Rl) Xl + {^Rl^Rl 

VR2) + b.c. , 



where 6 is the conjugate of related by = \T2^* , and 
a bare Majorana mass M is included. 



A. Charged lepton masses 

In our model, the 6*4 flavor symmetry is assumed to be 
spontaneously broken by the vacuum expectation values 
(VEVs) of Higgs scalars, i.e., {(pi) = w;, (6) = Ui, and 
(Xi) = ^i- One then arrives at the mass matrix of charged 
leptons as 

^ao - 2a2 63 62 \ 

Mp = \ 63 ao + %/3ai +03 61 (,10) 

&2 bi flo - V3ai +02/ 

where we have defined oo = a^VQ, (ai,a2) = (aiUi,aiU2), 
and bi = a2Ui (for i = 1, 2, 3). In general, all the param- 
eters in the mass matrix can be complex, while in case 
of CP-conservation, there are totally six real parameters 
in All. For simplicity, we will take all the parameters to 
be real, but comment later on the most general case with 
CP-violating effects. 

According to Eq. pop. the contributions from cpi 
merely affect the diagonal entries, whereas 6 appear in 
the off-diagonal elements. In the limit ai ^ bi, Mi ap- 
proximates to a nearly diagonal form, and the charged- 
Icpton masses are solely determined by a,; . Note that this 
is indeed a very realistic scenario if the VEVs of 6 are 
much smaller than those of (f>i. Explicitly, the sum of 
the VEVs has to be equal to the electroweak scale, i.e., 
Y,i |VEV,|^ ~ (174 GeV)^ Since cpo should also be re- 
sponsible for the generation of the top-quark mass, one 
may reasonably take vo ~ 174 GeV with all the other 
VEVs being much smaller than vq. In our model, we as- 
sume that vi,v2 ~ GeV and Ui ~ MeV. As a result, the 
eigenvalues of Mi are approximately given by oq — 2a2, 
ao 4- VSai + 02, and ao — VSai + a2, respectively. Com- 
pared to the charged-lepton masses, one immediately ob- 
tains 



1 / 

ao — 2 (me + m^ + mr) , 



ai 



2^3 



= (to^ - 77v) , 



a2 ~ - (m^ -t- rrir - 2me 



(11) 
(12) 

(13) 



In addition, the diagonalization matrix for M( is nearly 
an identity matrix, i.e., 14 ~ /. 



B. Neutrino mass matrix 

Since there are only two right-handed neutrinos in the 
MSM, the Dirac mass of neutrinos is a 3 x 2 matrix, viz. 



'2X1 
Md = I -X2 VIX2 



(14) 
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where Xi = for i = 1,2,3. The right-handed neu- 
trino mass matrix in our model is given by 



A + C 
B 



B 

A-C 



(15) 



where A = Mi, B ~ fiiXi, and C = /3ia;2, respectively. 
In case of Mo, Eq. ([T5|) leads to the masses of 

right-handed neutrinos as 



M- 



1,2 



A ± + C2 



(16) 



Note that, as aforementioned, the natural scale of AId 
relies on the VEVs implying Xi ^ ©(MeV). This in 
turn helps us to estimate that the right-handed neutrino 
masses should be around 0(10^) TeV, which turn out to 
be beyond the scope of forthcoming collider experiments. 
In case that certain fine-tuning is involved in the seesaw 
formula, e.g., the structural cancellation, one can, at least 
in principle, bring the masses of right-handed neutrinos 
down to the electrowcak scale, although the naturalness 
of such low-scale right-handed neutrinos seems question- 
able.^ 

By using the standard seesaw formula, i.e., TOj/ = 
—MdMj^^M]j, we obtain the light neutrino mass ma- 
trix as 

= mo X 

^262 - 2 (1 + ySei - £2)^1 (1 - \/3ei - £2)^2^ 

~ -(2 + x/3£i+£2)r2 (l + 2£2)rir2 1 (,17) 



[Viti - £2 - '2)rlj 



where 



Too 



2XlA 



(yl2 - 52 _ (^2) 



(18) 



and the parameters £ and r are defined by £1 = B/A, 
£2 = CI A, ri = X2/X1, and = X^/X^. 



C. Lepton flavor mixing 

The light neutrino mass matrix mj, is symmetric, and 
thus can be diagonalized by means of a unitary matrix Vn 
as Vjm^V* = diag(mi, 77i2, TO3). The lepton fiavor mix- 
ing matrix U which links the neutrino mass cigenstates 
with their fiavor eigenstates is then given by 



(19) 



^ Note that, the reaUzation of the TcV minimal seesaw model turns 
out to be more natural compared to the typical low-scale type-I 
seesaw model, since the light neutrino masses could be protected 
by certain underlying symmetries and hence do not suffer from 
large radiative corrections 1431. 



where the last approximation follows since we have taken 
the charged-lepton mass matrix to be nearly diagonal. In 
the standard (i.e., CKM-like) parametrization one has 

U = R23PsRi3Ps~^Ri2Pai , (20) 

where Rij correspond to the elementary rotations in the 
ij = 23; 13; and 12 planes (parametrized in what follows 
by three mixing angles c^^ = cos9^j and s^j = sm9^j), 
Ps = diag(l, l,e''^), and Pm = diag(e'"i/2^ ^1^2/2^ 
tain the Dirac and Majorana CP-violating phases, re- 
spectively. 

In order to get the explicit expression of U, a fully di- 
agonalization of rrii, is involved, and the results are rather 
tedious. However, since the to^ is of rank 2, there exists 
a eigenvector k = {rir2,r2,ri)'^ satisfying m^k = 0. In 
case that the light neutrino mass spectrum is inverted 
hierarchy (i.e., 7712 > mi 3> m^), k corresponds to the 
third column of U. Compared to Eq. ((20|) . we obtain 



tan 693 = — , 
tan 6^13 



(21) 
(22) 



In view of the experimentally measured maximal atmo- 
spheric angle and small reactor mixing angle, the relation 
ri ~ r2 ^ 1 has to be fulfilled. The two non-vanishing 
masses are then approximately given by 



TO2 



Too 



Too 



1 - £2 + v/(l - £2)2 +0{n,r2), (23) 



1 - £2 - - £2)2 +0{n,r2). (24) 



No matter what value of £2 one chooses, it is not pos- 
sible to let the two masses to be nearly degenerate (i.e., 
nil — TO2), which is indeed required for the inverted mass 
hierarchy case. Therefore, by analyzing the eigenvector 
of mi/, we can conclude that the inverted light neutrino 
mass hierarchy is not compatible with the model. 

Henceforth, we shall concentrate on the normal hier- 
archy case, namely mi < m2 <C ma. Here, we show one 
interesting example, in which the tri-bimaximal mixing 
pattern (i.e., 6I12 = 35.3°, 6*23 = 45° and 6*13 = 0) is 
predicted. Concretely, we make the assumptions that 
ri = r2 = 2 and £1 = 0. Equation (|25p now reduces to 

^2£2-2 2(1 -£2) 2(1 -£2) 
m„ = mo| - -4(2 + £2) 4(1 + 2£2) | . (25) 

-4(£2 + 2) 



One observes from Eq. ([25]) that, with the assumptions 
above, a /x — r symmetry appears in m^, which generally 
predicts a maximal atmospheric mixing angle, i.e., 623 ~ 
45°, and a vanishing ^13. It is then easy to prove that 
the diagonalization matrix of rrii, takes exactly the tri- 
bimaximal mixing form, i.e.. 



Utb 



\ 



(26) 
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while the light neutrino masses are given by 



The neutrino mixing angles are then modified to 



mi = , 

m2 = 6mo(l - 62) , 
ma = 12mo(l + £2) 



(27) 
(28) 
(29) 



Consequently, both the tri-bimaximal mixing and the 
normal neutrino mass spectrum (mi < m2 ma) are 
accommodated. 

Furthermore, if wc relax the assumptions on and e^, a 
deviation from the tri-bimaximal mixing can be achieved. 
Fox example, in the case ei ^ 0, the light neutrino mass 
matrix can be written as 



^0 

m„ = moC/TB I 6(e2 - 1) -6V2ei 



^0 - -12(1 + 62); 

UtbR23{0) diag(0,m2,m3) RI^Ut^^ (30) 



with 



m2 



ma 



3mo (^3 + 62 - ^(1 + 362)2 +8e2 

3mo (^3 + 62 + ^^(1 + 362)2 +862 



and 



sin 26* = 



2^/26l 



V/(1 + 362)2 + 86^ 



(31) 
(32) 

(33) 



S12 



S23 



Sl3 



^/i ^/3 
1 1 

sin 9 . 

V3 



sin^ - 
2 



(34) 
(35) 
(36) 



According to the above equations, the deviations of O^j 
from their exact tri-bimaximal values are correlated by 6*, 
and in the limit — )- (or effectively ei — )- 0), the exact 
tri-bimaximal mixing will be reproduced. Note that, the 
correction to S12 is proportional to sin^ |, and thus is 
strongly suppressed for a small 0. Therefore, ^12 is rather 
stable against ei corrections [48ll49{. 



D. Higgs potential 

Now that the previous discussions rely on the VEVs 
of the scalar fields, we are coming to the question of the 
possible Higgs potential and its minima. Apart from the 
SU (2) singlets Xi our Higgs setup is essentially the same 
as the Higgs sector considered in Ref. [l^l, where only 
SU{2) doublets are introduced. We thereby only show 
the Higgs potential parts involving Xi, viz.. 



14 



-A^i {xl + xl) +/"2 (3x?X2 - xl) +^^1 (x? +^2)^ +^2 [(xiX2 +X2X1) + (xl - xl)]'^ 



Pi 
Pi 



'^o'^o (xi - X2) +^2^0^X1 +0^X2 j + Pa (0501X1 + 0^02X2 + h.c.^ 
0S01 (XiX2 + X2X1) + 0S02 {xl - xl) + h.c. 



402) (X?+X2) +£2 



0l02 + 020l) (XlX2 + X2X1) + (0l01 - 0202) (Xl - X2) 



£3 



0l02 + 020l) Xl + (0l01 - 0202) X2 

2 . t 2\ 

01X1 - 02X2 I 



£4 



^'IXl +02X2 



£5 



4x2 + 4>\xi 



4X2 + 02X1 



ki (ela + ek2 + eka) {xl + X2) + ^2 [Vi (ek2 - eks) (xiX2 + X2xi) + (e2k2 + eks - seki 



{xl xl) 



^3(41^X1!^+ %/36xi+6X2 + V3CaXi-6 



^3X2 



fc4(4|CiXi|^+ ^3^X2- 6X1 + ^/3C3X2+6Xl 



^ (C2k2 - eks) Xl + (-2eki + ek2 + eks) X2 



(37) 



r 



Compared to the Higgs potential in Ref. [20| , there are in VEV structure described in the previous analysis can be 
total 21 more parameters. So we are confident to arrive easily satisfied. Furthermore, we did not discuss in de- 
al the suitable minima of the Higgs potential, and the tail the Higgs spectrum, which may result in the flavor 
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FIG. 1: The allowed parameter regions in the ri — r2 plane 
(upper plot) and the ei — £2 plane (lower plot). 

changing neutral currents as well as lepton flavor violat- 
ing problems. However, such problems commonly occur 
in models with more than one Higgs doublets, and can be 
ignored if the flavor changing Higgs arc all heavier than 
a few TeV. 



III. NUMERICAL ILLUSTRATIONS 

We proceed to the numerical illustrations. The in- 
put values for the neutrino parameters are taken from 
Ref. [soj . For example, in the normal hierarchy case, the 
mass-squared differences measured in atmospheric and 
solar neutrino experiments read 

Aml^ = (7.12 - 8.13) x 10"^ eV^ , (38) 
Am^i = (2.18 - 2.73) x 10"^ eV^ , (39) 

while the allowed ranges of three mixing angles are 

sin^^iz = 0.27 ~ 0.37, (40) 
sin2 6'23 = 0.39 - 0.64, (41) 
sin2 6'i3 < 0.04, (42) 



7 r 
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2 I , , , , , 1 

-2-1012 

El 

FIG. 2: The allowed regions of the ratio M\/M2 with respect 
to ei. 

at 3cr confidence level. Note that, there are slightly dif- 
ferences between the fitted parameters in the inverted 
and normal hierarchies. In the normal hierarchy case, 
the above mass-squared differences correspond to the al- 
lowed range of the mass ratio 5.4 < 7713/7/12 < 5.8. 

In our numerical analysis, we do not make any assump- 
tions on the model parameters, and randomly choose the 
values of r^, and ttiq. The predicted neutrino mixing 
angles and masses (in the normal hierarchy case) are then 
compared with Eqs. ([55|) - P^ . while the allowed param- 
eter spaces of and are shown in Fig. [TJ From the 
upper plot, one observes that the allowed regions of ri 
and r2 are symmetric, which is actually resulted from the 
V — T symmetry in the neutrino mass matrix. In addi- 
tion, none of or r-i can be zero, while ri ~ r2 ~ 2 is 
quite favored according to the numerical results. In the 
lower plot; ei = is allowed but £2 = is not, indicating 
that X is required in order to fit the experimental data. 
Furthermore, for a fixed value of ei, there are two al- 
lowed regions for £2 corresponding to £2 > 1 and £2 < 1, 
respectively. 

Since the right-handed neutrino masses are also corre- 
lated to Ti , we present in Fig. [2] the predicted mass ratio 
between two right-handed neutrinos. One reads from the 
figure that the mass ratio is generally larger than 2 show- 
ing that the resonant leptogensis mechanism[5lj may not 
simply apply to this model. ^ 

Now we turn to the special case with the assumption 
ri = r2 = 2, namely, a /i — t symmetry exists in the 
neutrino sector. The allowed parameter regions of dij 
and M1/M2 arc illustrated in Fig. [3] As we expected, 
there exist strong correlations between three mixing an- 
gles according to the upper and middle plots. This is 



^ The right-handed neutrinos are degenerate in Ref. [20l | since their 
masses are originated from a bare Majorana mass term, whereas 
in our model, due to contributions from Xi a mass spUtting be- 
tween Ml and M2 is included. 



6 



60 I , , , , , 1 

50 

40 - 

°^ 30 - 

CD 

20 
10 - 

I < < < < < 

1 2 3 4 5 6 




°^ 30 - 

CD 

20 
10 

I ' ' ' ' ' 

1 2 3 4 5 6 

3.5 I < < < < < 1 

3 - 

2.5 

« 2 - 

1.5 - 

1 - 
0.5 - 

I ' ' ' ' ' 

1 2 3 4 5 6 

FIG. 3: The allowed parameter regions of 6ij and M\/M2. 

in good agreement with our analytical results aforemen- 
tioned since the mixing angles are connected by a single 
parameter 6. The most severe constraint comes from 023, 
and its experimental allowed range can be fulfilled. As 
for 013, an upper bound ^13 < 5° can be obtained. As has 
been shown, 612 is confined to it's tri-bimaximal mixing 
value, and rather stable compared to the two other mix- 
ing angles. In the particularly interesting limit ei = 0, 



0.01 




2 4 6 8 10 12 
FIG. 4: The allowed parameter regions of \mee\ and 613. 



the exact tri-bimaximal mixing pattern will be repro- 
duced. Finally, from the lower plot, we also find that 
the right-handed neutrino mass spectrum should be hi- 
erarchical, e.g., M1/M2 ~ 3. 

We stress that our discussions are based on the as- 
sumption of real Yukawa couplings as well as scalar 
VEVs, whereas in the most general situation, both of 
them could be complex. In the presence of CP-violating 
effects, the imaginary parts of the model parameters 
A, B, C and X could significantly change the predic- 
tions addressed here, and we therefore study it further. 
Since the neutrinoless double beta decay process rely on 
the Majorana feature of light neutrinos, we illustrate 
in Fig. |4] the allowed ranges of the effective mass, i.e., 
™ee = ''^iUei-i ^-^d 013 in the presence of CP violation. 
The model parameters are the same as those in Fig. [31 
except that we allow them to be complex. One observes 
from the plot that any value of 0i3 satisfying the cur- 
rent experimental constraint can be achieved, whereas 
there exist strong constraints on |mee|. in particular for 
a smaller 0i3, indicating potentially attractive signatures 
in future non-oscillation experiments. 



IV. CONCLUSION 

In this work, we presented a minimal seesaw model 
based on the discrete iS'4 flavor symmetry. In our model, 
besides the SM fermion content, two right-handed neu- 
trinos are introduced transforming as an S'4 doublet. The 
structure of the model is minimal in the sense that there 
are at most two massive light neutrinos which are in- 
deed required to account for the observed solar and at- 
mospheric neutrino oscillations. The number of model 
parameters are reduced greatly compared to the ordinary 
type-I seesaw, and thus allow us to make useful predic- 
tions on the neutrino parameters. After carefully ex- 
ploring the parameter spaces, we found that the inverted 
neutrino mass hierarchy is ruled out whereas the nor- 
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mal hierarchy can be well accommodated in this frame- 
work. In particular, the tri-bimaximal mixing pattern 
can be naturally obtained from simple assumptions on 
the model parameters, while the deviation of three mix- 
ing angles from their exact tri-bimaximal mixing values 
are correlated by a single model parameter. In addition, 
the right-handed neutrinos feature a hierarchical mass 
spectrum, i.e., the ratio between right-handed neutrino 
masses is generally larger than 2. 

Note that, in the current discussions, we have ignored 
the CP- violating effects, since there is yet no direct ex- 
perimental information on leptonic CP violation. How- 
ever, in the most general case, the CP-violating phases 
can be easily included since all the coefficients of Yukawa 
couplings as well as the VEVs could in principle be com- 
plex. In fact, the CP-violating effects are very crucial in 
order to explain the baryon asymmetry of the Universe 
via thermal leptogenesis mechanism in the seesaw mod- 
els [S^l ■ In addition, a dirac CP-violating phase may also 
be searched for at future long-baseline neutrino oscilla- 
tion experiments. 

Finally, we stress that the right-handed neutrinos may 
not be necessarily heavy, e.g., their masses could be lo- 
cated around keV scales. One may wonder that, in the 
mass range Mi ~ keV (i.e., the right-handed neutrinos 



are sterilized), if the right-handed neutrinos could be 
viewed as warm dark matter so as to explain simultane- 
ously the neutrino mass generation and the dark matter 
puzzle. Unfortunately, this is not possible in the current 
model, since the stability of keV right-handed neutrinos 
on the cosmic time scale requires the mixing between 
sterile and active neutrinos to be smaller than 10~^, 
which leads the mass scale of light neutrinos to be about 
10~^ eV [S^l- Such tiny neutrino masses are in conflict 
with neutrino oscillation experiments. Possible variations 
extending the MSM may provide successful warm dark 
matter candidate, (e.g., an additional light right-handed 
neutrino transforming as a singlet under 6*4), which are 
however beyond the scope of current work. 
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